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MERLIN imaging of the maser flare in Markarian 348 
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ABSTRACT 

MERLIN images of Mrk 348 at 22 GHz show H 2 maser emission at 0.02 - 0.11 Jy, 
within ~ 0.8 pc of the nucleus. This is the first direct confirmation that molecular 
material exists close to the Seyfert 2 nucleus. Mrk 348 was observed in 2000 May 
one month after Falcke et al. (2000) first identified the maser in single-dish spectra. 
The peak maser flux density has increased about threefold. The masing region is 
^ 0.6 pc in radius. The flux density of radio continuum emission from the core has 
been rising for about 2 years. The maser-core separation is barely resolved but at 
the 3cr significance level they are not coincident along the line of sight. The masers 
lie in the direction of the northern radio lobes and probably emanate from material 
shocked by a jet with velocity close to c. The correlation between the radio continuum 
increase and maser flare is explained as arising from high level nuclear activity through 
a common excitation mechanism although direct maser amplification of the core by 
masers tracing a Keplerian disc is not completely ruled out. 

Key words: masers - galaxies: individual (Markarian 348) - galaxies: nuclei - galaxies: 
Seyfert - radio lines: galaxies - radio continuum: galaxies 
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1 INTRODUCTION 

Water maser emission has been detected from galaxies up 
to ~ 100 Mpc from Earth. These supermasers are hundreds 
of times more luminous than the brightest masers in star- 
forming regions in our own Galaxy (Cohen 1998). Such 
bright spectral lines are ideal targets for very long baseline 
radio interferometry with micro-arcsec accuracy and with 
Si 1 km s _1 velocity resolution. This is the only way to 
image distant galaxies directly on sub-pc scales. 

H2O supermasers are found exclusively in about 5% of 
type 2 Seyfert and LINER galaxies (Braatz et al. 1997). 
The unified scheme predicts that these active galactic nuclei 
(AGN) are obscured by a molecular torus viewed edge on. 
This provides high column densities for maser amplification 
(Kartje et al. 1999). Such masers have been used to measure 
the parameters of nuclear discs in Keplerian rotation (e.g. 
Miyoshi et al. 1995). It is less obvious why the selection 
effect applies to supermasers which appear to be associated 
with a jet instead of a circumnuclear disc (e.g. Claussen 
et al. 1998). These masers could originate from interstellar 
material shocked by the jet (Elitzur et al. 1992). 

Falcke et al. (2000) reported the discovery of a very lu- 
minous H2O maser in Mrk 348 during a radio flare of the 
AGN. Mrk 348 is a well studied Seyfert 2 at a redshift of 
0.015 (Huchra et al. 1999), with broad emission lines in 
polarized light (Miller & Goodrich 1990). However its host 
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is an SO galaxy at an angle of inclination of only « 16° 
(Simkin et al. 1987). Ground-based (Simpson et al. 1996) 
and HST (Falcke et al. 1998) imaging have revealed a dust 
lane crossing the nucleus. It has a high x-ray-absorbing col- 
umn depth of JVh = 10 27 ' 1 m -2 towards the nucleus (War- 
wick et al. 1989). These observations suggest the presence 
of an obscuring torus in Mrk 348, but no molecular or HI 
absorption has been detected so far (Gallimore et al. 1999; 
Falcke et al. 2000). 

In this paper all velocities are given relative to the local 
standard of rest (Vlsr) in the radio convention. The sys- 
temic Vlsr of Mrk 348 is 4435 < V sys < 4480 km s" 1 from 
Hi emission line observations (Bottinelli et al. 1990) and 
imaging (Simkin et al. 1987). We adopt Ho = 75 km s -1 
Mpc -1 , so Mrk 348 is at a distance of ~ 60 Mpc, where 1 
mas = 0.29 pc. 

What makes Mrk 348 stand out among Seyfert galaxies 
is its bright and variable inverted-spectrum radio nucleus 
(Unger et al. 1984). Jets at position angles (p. a.) 170° 
and 30° were observed with the EVN at 1.4 GHz (Neff 
& de Bruyn A G 1983) and MERLIN at 5 GHz (Unger 
et al. 1984). Halkides et al. (1997) first resolved the cen- 
tral part of Mrk 348, using the VLBA at 15 GHz, into two 
components separated by ~ 0.3 pc at p. a. ~ 90° to the 
larger 1.4 GHz jet. Ulvestad et al. (1999) measured the sub- 
relativistic separation rate of these two components at the 
same frequency. They also noted a rise in continuum flux 
from 120 to 570 mjy between 1997.10 and 1998.75. Using 
MERLIN we can detect mjy radio continuum with a surface 
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brightness <; 5 x 10 4 K and locate bright masers with sub-pc 
relative positional accuracy. We use this to investigate the 
relationship between the maser and continuum flares and 
their origins in the core or in a nascent jet. 



2 OBSERVATIONS AND DATA REDUCTION. 

We observed Mrk 348 on 2000 May 2 using MERLIN, which 
has a maximum baseline of 217 km, giving a beam size of 
12 mas. In order to obtain images of Mrk 348 as rapidly as 
possible we observed for a single 17 hr run in the maximum 
16 MHz bandwidth, which corresponds to a total velocity 
width of ~ 200 km s _1 and only covers the red-shifted half 
of the line seen by Falcke et al. (2000) plus sufficient line-free 
channels for continuum subtraction. We observed Mrk 348 
at a fixed frequency of v — 21891.6 MHz alternately with 
phase reference source J0057+3021 for 4 and 2 min respec- 
tively. 3C273, which had a flux of 21.09 Jy at that time 
(Terasranta, private communication), was observed for 40 
min and used to set the flux scale for all sources. 

Further data processing was performed using MPS 
(Greisen 1994). We applied the instrumental corrections 
from the calibrators and the phase reference source solu- 
tions to Mrk 348. We then adjusted the data to a constant 
velocity, putting Vlsr = 4641.6 km s _1 in channel 30 of 60 
usable channels, with a separation of 3.37 km s _1 . We av- 
eraged all data for Mrk 348 and made a CLEAN map using 
natural weighting of the visibility data, which gave a beam 
size of 31 mas x 21 mas. 

The three clean components (cc) above 3cr rm s from 
the initial map of Mrk 348 lay in a region 6 by 18 mas elon- 
gated north-south. We used these as a model for phase self- 
calibration, applied the solutions and examined the spec- 
trum of Mrk 348. The complex visibilities were vector aver- 
aged in time channel by channel. The flux density increased 
noticably with increasing frequency if we used the more 
southerly CC positions as the origin of phase. Fig. [j] shows 
the spectrum on the baselines to the Cambridge antenna us- 
ing the mean position of the maser emission (Table [l]) as the 
phase origin. The low-frequency end of the band appeared 
to be continuum only. We used the first 4.25 MHz of data for 
phase and amplitude self-calibration and applied the same 
corrections to the line-only data as we did for the continuum 
data. 

The multi-channel Mrk 348 data were averaged over 
every 0.75 MHz and Fourier transformed to make a 20- 
channel total emission dirty data cube. The average of the 
first 4.50 MHz was subtracted pixel by pixel from the whole 
cube to leave a line-only dirty data cube which was CLEANed 
to produce a CLEAN line-only cube. The line-only dirty data 
cube was subtracted from the total emission dirty data cube 
and the result was CLEANed to give a continuum-only cube. 
We also made a CLEAN total emission cube. 

We fitted 2D Gaussian components to emission above 
3<7 rms in each channel of every CLEAN cube to determine 
the position and peak flux S p . We also measured the posi- 
tion and flux density of components fitted to the total flux 
from Mrk 348 mapped using 15-MHz band-width prior to 
self-calibration and the 4.25-MHz bandwidth self-calibrated 
continuum map. The brightest continuum components was 
resolved after deconvolving the beam, so we could measure 



Frequency (GHz) 

21.89 21.888 21.886 
1 ' 1 ' 1 




. I I i i i i I i i i i I i i i i I i i i 

°4550 4600 4650 4700 

^lsr ( racl ' convention, km s~ 1 ) 

Figure 1. The total flux density of Mrk 348 in each channel. 
The arrows mark the frequencies of peaks observed by Falcke et 
al. (2000). 

its FWHM s and the integrated flux Si. The position and 
component size uncertainties (cr P os, er s ) are proportional to 
the beam-width divided by the signal-to-noise ratio allow- 
ing for the sparse baseline coverage (Condon et al. 1998; 
Richards et al. 1999). 

All maps are presented in total intensity. The flux scale 
should have < 10% error but MERLIN only has 5 antennas 
operating at 22 GHz, and the sparse uv coverage means it is 
possible for extended continuum emission to appear brighter 
and more compact than it actually is, although the peak po- 
sitions should be accurate. Moreover as only one side of the 
bandpass is line-free, errors in baseline subtraction or band- 
pass calibration may affect the maser flux measurements. 
The absolute position accuracy is ~ 50 mas, mostly due to 
uncertainty in the position of J0057+3021 (Wilkinson et al. 
1998). 



3 RESULTS 

A single patch of maser emission was detected in 11 
channel maps of the line-only data cube from 4543.8 to 
4655.0 km s _1 , shown by the contours in Fig. ^. These 
are in linear multiples of 3cr rms . The grey-scale shows the 
continuum-only emission above 3(j rms . The continuum peaks 
are marked N and S and the white crosses show their posi- 
tions, which agree in each channel to within 0.1 mas. 

The positions and flux densities of the peaks are given 
in Table |l| The peak of the maser emission is consistently 
offset from the continuum peak; its mean position M is 
2.7 ±0.7 mas north of S. The total angular size of the maser 
region is < 3 mas and the individual components are unre- 
solved. There is no significant systematic angular separation- 
velocity gradient greater than 2 mas in 111 km s _1 . S ap- 
pears to be resolved, and was fitted with a component of 
2.0 ± 0.2 mas FWHM and an integrated flux of 814 mjy 
giving a brightness temperature of 551 - 824xl0 6 K. N is at 
a position angle of —11° ± 2° with respect to S. 

The maps of continuum-only channels from 21884.35 to 
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Table 1. The positions and peak intensities S p of the 22 GHz continuum and maser peaks in Mrk 348 N, S and M. 

Peak R.A. Dec. o" P os (mas) S p (mjy beam -1 ) <r rms (mjy beam - x ) 

S 00 h 48 m 47!14575 +31° 57' 25'.'1128 0.2 744 6 

N 00 h 48 m 47!14552 +31° 57' 25^'1470 0.4 239 6 

M 00 h 48 m 47!14580 +31° 57' 25^1155 0.7 24 - 107 8 
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Figure 2. H2O maser emission from Mrk 348 is shown by the 
contours, at (-1, 1, 2, 3, 4...)x20 mjy beam . The grey scale 
shows continuum emission from to 1 Jy beam -1 . The circular 
restoring beam of 12 mas FWHM is shown in the top left panel. 
The white crosses show the continuum peaks N and S, size pro- 
portional to the peak flux. The axes are labelled in mas offset 
from S. Each panel is labelled with the Vlsr (radio convention). 



21888.60 MHz, the continuum+line channels from 21890.85 
to 21899.10 MHz and line-only channels from 21890.85 to 
21899.10 MHz were averaged from their respective data 
cubes. Flux measurements along slices through these aver- 
aged maps, at constant R.A., 2 mas wide, intersecting the 
position of S, are plotted in Fig. || This shows the offset of 
the maser peak from the brighter continuum peak. 



4 DISCUSSION 

These observations are the first direct detection of molecu- 
lar material in the nuclear region of Mrk 348. The existence 
of H2O implies shielding by a dusty medium with a col- 
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Figure 3. Flux density as a function of position along slices 
through the maser and continuum peaks. The dot-dashed line 
shows the slice taken from the average of continuum-only chan- 
nel maps. The dotted line shows the average of channel maps 
containing maser and continuum emisson. The solid line shows 
the averaged maser emission in these channels after continuum 
subtraction, note the different flux scale. 



umn density 10 26 — 10 27 m 2 and the conditions required 
for population inversion of the maser include a gas number 



densitity 
~ 10" 5 - 



10 1 



10 m , a fractional abundance of H2O 



10 4 and a temperature > 250 K (Kartje et al. 



1999). Conditions can be more tightly constrained depend- 
ing on the association of the masers with a disc or a jet. 

We assume that the brightest 22-GHz continuum com- 
ponent S contains the core and 2-mas northern jet at posi- 
tion angle 15° found by Ulvestad et al. (1999) at 15 GHz. 
N is in a similar direction to the more extended jet detected 
at 1.4 and 5 GHz (Neff & de Bruyn A G 1983; Unger et al. 
1984). 

We consider three possible models for the relationship 
between the maser flare and the continuum flare: 

(i) The masers trace a small warped disc and the maser 
emission follows a Keplerian rotation law, and the masers 
directly amplify the continuum emission. 

(ii) The masers are unsaturated and lie in a symmetric 
Keplerian nuclear disc. Perturbation of the disc creates spi- 
ral shocks (Maoz & McKee 1998) that are also responsible 
for the continuum flare. 

(iii) The radio continuum flare is associated with the ejec- 
tion of material in the direction of the northern jet. The 
masers arise from the ISM where it is shocked by the jet. 

Below we expand on the models and explain why we prefer 
or reject each. 

The ratios of the peak fluxes ApribMay are 27:77 at 
4604 km s" 1 and 26:107 at 4544 km s" 1 , a 3 - 4-fold increase. 
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Although further observations are essential in order to pin 
down the source timescale, such rapid variability may sug- 
gest that the individual masing clouds are less than 1 light- 
month (0.025 pc) in diameter. This would then be about half 
the size of the 15-GHz core (Ulvestad et al. 1999) and so 
the maser emission that is beamed in our direction is com- 
pact enough to originate from maser clouds that overlap our 
line of sight to the core. However the two arrows in Fig. fil 
indicate the velocities of the peaks that Falcke et al. (2000) 
observed and which are close to the peaks in our data, but 
~ 100 km s^ 1 more redshifted than V sys . Falcke et al. (2000) 
did not detect any maser emission at Vlsr < U sys . So ex- 
planation (1) is inconsistent with a Keplerian disc unless 
the masing material is infalling. Moreover the maser peak 
is consistently misaligned with S in every channel in Fig. [| 
The 3cr position error boxes for M and S (Table [j]) are too 
close to rule out direct maser amplification of the continuum 
peak, but we consider other geometries are more probable. 

Maoz & McKee (1998) found a tendency for the red- 
shifted emission from disc supermasers around AGN to be 
brighter and explained this using a model of spiral shocks 
within a Keplerian disc. This mechanism requires pertur- 
bation of the disc, which could be connected to the event 
which caused the continuum flare. Such a disc should be 
symmetric about the nucleus with a well-defined position- 
velocity gradient. If it is nearly edge-on the rotation velocity 
is Vmax - V sys « 200 km s" 1 . The MERLIN results show the 
velocity gradient is < 2 mas in 111 km sj -1 corresponding 
to a disc radius 1.1 pc. M is 0.8 ± 0.2 pc north of S. This 
is consistent with the 2 yr time-lag between detection of 
the maser and continuum flares which implies a separation 
of <; 0.6 pc (Falcke et al. 2000). This implies that if the 
masers lie in a symmetric disc it is elongated north-south, 
parallel to the radio jets. This is unlikely and so model (2) 
is also improbable. 

Model (3), in which the masers originate in a shock pro- 
duced by the northern jet, is possibly due to the misalign- 
ment of the nuclear disc with respect to the host galaxy. 
If this material was ejected from the core when the contin- 
uum flare commenced 2 years prior to the maser flare 
this implies speeds near c . The only previous detection of 
a relativistc Seyfert jet is in III Zw 2 (Brunthaler et al. 
2000). However Ulvestad et al. (1999) measured a jet speed 
of ~ 0.07c for Mrk 348 on similar scales suggesting the jet 
power is rapidly dissipated once it reaches the shocked re- 
gion. 

If the maser flare is not simply direclty amplifying the 
radio continuum flare this indicates that there is some com- 
mon excitation effect, possibly some sort of high level nuclear 
activity. If the masers the masers are found along the radio 
jets, as is the case for NGC 1052 (Claussen et al. 1998), this 
mechanism for the correlation between the evolution of the 
maser flare and the radio flare is an important tool to study 
the jet-ISM interactions. 

We have observed the newly discovered megamaser 
Mrk 348 with MERLIN at 22 GHz in order to study the 
masers at higher resolution, test whether there is a corre- 
lation between the continuum and maser flux density and 
investigate their origin in the core or jet, which could give 
us some more insight to the maser excitation mechanism. 
Further imaging observations of the whole line will distin- 



guish between these possibilities and elucidate the origins of 
the flare. 
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